BOUNDARY DRIVEN KAWASAKI PROCESS WITH LONG 
RANGE INTERACTION: DYNAMICAL LARGE DEVIATIONS 
AND STEADY STATES 



MUSTAPHA MOURRAGUI AND ENZA ORLANDI 

Abstract. A particle system with a single locally-conserved field (density) in 
a bounded interval with different densities maintained at the two endpoints of 
the interval is under study here. The particles interact in the bulk through a 
long range potential parametrized by f) > and evolve according to an exclu- 
sion rule. It is shown that the empirical particle density under the diffusive 
scaling solves a quasi-linear integro-differcntial evolution equation with Dirich- 
let boundary conditions. The associated dynamical large deviation principle 
is proved. Furthermore, for /3 small enough, it is also demonstrated that the 
empirical particle density obeys a law of large numbers with respect to the 
stationary measures (hydrostatic). The macroscopic particle density solves a 
non local, stationary, transport equation. 



1. Introduction 

Over the last few years there have been several papers devoted to understand- 
ing the macroscopic properties of systems out of equilibrium. Typical examples 
are systems in contact with two thermostats at different temperature or with two 
reservoirs at different densities. A mathematical model is provided by interacting 
particles performing local reversible dynamics (for example reversible hopping dy- 
namics) in the interior of a domain and some external mechanism of creation and 
annihilation of particles at the boundary of the domain, modeling the reservoirs. 
The full process is non reversible. The stationary non equilibrium states are char- 
acterized by a flow of mass through the system and long range correlations are 
present. We refer to [3] and [6] for two reviews on this topic. We study a model 
with such features but in a regime where phase separation might occur at equilib- 
rium for the underlying reversible model. To this aim we consider in the interior of 
the domain a reversible dynamics (Kawasaki dynamics) constructed trough a mean 
field interaction (Kac potential), see below for more details. This is the first time, 
to our knowledge, where both, long range dynamics in the bulk and creation and 
annihilation of particles at the boundary of the domain, are taken into account. 

The particle models we consider are dynamic versions of lattice gases with long- 
range Kac potentials, i.e. the interaction energy between two particles, say one at 
x and the other at y (x and y are both in Z d ), is given by Jn(x, y) = N~ d J(jf, 
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where J is a smooth symmetric probability kernel with compact support and ./V 
is a positive integer which is sent to oo. The equilibrium states for these mod- 
els have been investigated thoroughly [12], [18] and [22], and have provided great 
mathematical insight into the static aspects of phase transition phenomena. The 
dynamical version of these lattice gases in domain with periodic boundary condi- 
tions (reversible dynamics) has been analyzed in [8], [9], [10] and [1]. This paper 
starts by studying the dynamics of these systems in a bounded interval with reser- 
voirs (non-reversible dynamics). We investigate their qualitative behavior in the 
range of the parameter when at equilibrium phase transition is present. Let us 
describe informally the dynamics. We consider a one dimensional lattice gas with 
particle reservoirs at the endpoints. The restriction on the dimension is done only 
for simplicity. Given an integer N > 1 at any given time each site of the discret 
set {— N, . . . , N} is either occupied by one particle or empty. The interaction en- 
ergy among particles is given by a modified version of the Kac potential Jn and 
it is tuned by a positive parameter /?. The modification of the Kac potential, see 
(2.1), takes into account that the particles are confined in a bounded domain. In 
the bulk each particle jumps at random times on the right or on the left nearest 
neighborhood, if the chosen site is empty, at a rate which depends on the particle 
configuration through the Kac potential. When (3 = we have the simple exclusion 
process. At the boundary sites, ±iV, particles are created and removed for the local 
density to be < p- < p+ < 1. At rate p± a particles is created at ±N if the site 
is empty and at rate 1 — p± the particle at ±iV is removed if the site is occupied. 

The dynamics described above defines an irreducible Markov jump process on a 
finite state space, i.e. there is a strictly positive probability to go from any state 
to another. By the general theory on Markov process [17], the invariant measure 
/j s i"' is unique and encodes the long time behavior of the system. Let P ,f„, be 

the stationary process. Since p s i at is invariant, the measure V^'JtLt is invariant with 

respect to time shifts. The measure P J^ at is invariant under time reversal if and 

only if the measure p s ^ at is reversible for the process, i.e if the generator of the 
process satisfies the detailed balance condition with respect to p s ^ at . In the case at 
hand p s ^ at is stationary and reversible only if (3 = and p + = p~ = p. In such case 
the invariant, reversible measure is the Bernoulli product measure with marginals 
p. When f3 > and /or p + ^ p~ then p s ^ at is not reversible. In such case the 
corresponding process is denoted non reversible. We shall consider the latter case. 
The lattice space is scaled by and the time by N 2 (diffusive limit) while the 
behavior is studied, as N f oo, of the empirical density of the particles evolving 
according to the dynamics described above. 

To prove the hydrodynamic behavior of the system, we follow the entropy method 
introduced by Guo, Papanicolaou, and Varadhan [11]. It relics on an estimate on 
the entropy of the state of the process with respect to a reference invariant state. 
The main problem in the model considered here is that the stationary state is not 
explicitly known. We have therefore to consider the entropy of the state of the 
process with respect to a state which is not invariant, for instance, with respect to 
a product measure with a slowly varying profile. Since this measure is not invariant, 
the entropy does not decrease in time and we need to estimate the rate at which it 
increases. These estimates on the entropy are deduced in Subsection 3.1 
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It results that the local empirical particle density, in the diffusive limit, is the 
solution of a boundary value problem for a quasilinear integro differential para- 
bolic equation, see (2.9). In addition to this, it is demonstrated that when /3 is 
small enough, then the empirical particle density obeys a law of large numbers 
with respect to the stationary measures (hydrostatic). The result is obtained char- 
acterizing the support of any limit points of P . This intermediate result holds 
for any /? > 0. Then it is shown that with /3 small enough, the stationary solution 
of the boundary problem (2.9) is unique, it is a global attractor for the macroscopic 
evolution with a decay rate uniform with respect to the initial datum. This holds 
only for j3 < (3q where (3$ > depends on the diameter of the domain and on the 
chosen interaction J. Namely the quasilinear non local parabolic equation does not 
satisfy a comparison principle, which is the main tool used in previous papers, see 
for example [20] and [7], to show the hydrostatic. For value larger than (3q we are 
not able to show the uniqueness of the stationary solution of the boundary value 
problem (2.9). We stress that /?o < /3 C where f3 c is the value above which phase 
segregation occurs at equilibrium; with the choice made of the parameters j3 c = j , 
see page 1712 of [9]. Further, we prove the dynamical large deviations for the em- 
pirical particle density. The large deviation functional is not convex as function of 
the density, it is lower semicontinuous and has compact level sets. Since the large 
deviation functional is not convex and the underlying dynamics does not satisfy any 
comparison principle, care needs to be taken to prove the lower bound. The basic 
strategy to show the lower bound consists in obtaining this bound for smooth paths 
and then applying a density argument. The argument goes as follows: Given a path 
p with finite rate functional I{p) one constructs a sequence of smooth paths p n so 
that p n — > p in a suitable topology and I{p n ) — > I(p)- When the large deviation 
functional is convex, this argument is easily implemented. In our case, because of 
the lack of convexity we modify the definition of the rate functional declaring it 
infinite if a suitable energy estimate docs not hold. In this way the modified rate 
functional when finite provides the necessary compactness to close the argument. 
This method has been developed in [21] and we adapted to our model. The mod- 
ification of the rate functional helps in showing the lower bound but makes more 
difficult the upper bound. One needs to show that the energy estimate holds with 
probability super exponentially close to one. A similar strategy was applied in [20] 
and [4]. 

In a recent paper, De Masi et alii, [5], constructed, in the phase transition regime, 
stationary solutions of a boundary value problem equivalent to (2.9) in which the 
density p is replaced by the magnetization m = 2p — 1. They did not study the 
derivation of the boundary value problem from the particle system and they did 
not inquire about the uniqueness of stationary solutions. They investigated the 
qualitative behavior of constructed stationary solutions of (2.9) as the diameter of 
the domain goes to infinity. They proved, for the constructed solution, the validity 
of the Fourier law in the thermodynamic limit showing that, in the phase transition 
regime, the limit equilibrium profile has a discontinuity (which defines the position 
of the interface) and satisfies a stationary free boundary Stefan problem. 

The paper is organized as follows: The precise feature of the model, notations 
and results are stated in Section 2, In Section 3, some basic estimates needed along 
the paper are collected. In Section 4 the hydrodynamic and the hydrostatic limits 
are shown. Section 5 is split into 5 subsections and deals with dynamical large 
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deviations. We prove in Section 6 the weak uniqueness of the solution of the quasi- 
linear integro-differential equation. Furthermore, for /3 small enough, it is shown 
that its stationary solution is unique and it is a global attractor in L 2 . 

2. Notation and Results 

Fix an integer N > 1. Call A N = {-N,...,N} and T N = {—N, +N} the 
boundary points. The sites of Ajy are denoted by x, y and z. The configuration 
space is S N = {0, 1} An which we equip with the product topology. Elements of Sn, 
called configurations, are denoted r\ so that r)(x) <E {0, 1} stands for the number of 
particles at site x of the configuration r\. 

We denote A = (—1,1) (A = [—1,1]) the macroscopic open (closed) interval, 
r = { — 1, 1} its boundary and u G [—1, 1] the macroscopic space coordinate. 

2.1. The interaction. To define the interaction between particles we introduce 
a smooth, symmetric, translational invariant probability kernel of range 1, i.e 
J(u, v) = J(0, v — u) for all u, v £ R, J(0, u) = 0, for all |u| > 1, and / J(u, v)dv = 1, 
for all u £ R, 1 . When (u, v) £ A x A we define the interaction J neum (u,v) im- 
posing a reflection rule: u interacts with v and with the reflected points of v where 
reflections are the ones with respect to the left and right boundary of A. For this 
reason it is referred to as "Neuman" interaction. More precisely, we define for u 
and v in A 

J neum (u, v) := J(u, v) + J(u, 2-v) + J(u, -2 - v), (2.1) 
where 2 — v is the image of v under reflections on the right boundary {1} and —2 — v 
is the image of v under reflections on the left boundary { — 1}. By the assumption 
on J, J neum (u,v) = J neum (v,u) and / J neum {u, v)dv = 1 for all u e A, see Lemma 

3.1. We defined the interaction (2.1) by boundary reflections only for convenience. 
It has the advantage to keep J neurn a symmetric probability kernel. This choice of 
the potential has been done already in [5] . 

The pair interaction between x and y in Ajv is given by 

J N (x,y)=N- 1 J~(^,^). 
The total interaction energy among particles is given by the following Hamiltonian 

H n(v) = -\ J2 MwWxWv) ■ (2-2) 

2.2. The dynamics. We denote by rj x - y the configuration obtained from 77 by 
interchanging the values at x and y: 

{v(y) if z = x 
V (x) if z = y 
r](z) if z^x,y, 

and by a x rj the configuration obtained from 77 by flipping the occupation number 
at site x: 

{(jXrj){z):= ^-r,{x) iiz = x 
I r](z) if z/i. 



^Onc could take the interaction J so that for all u 6 R, J J(u,v)dv = a with a > 0. The 
only difference, with the case at hand, is that the the underlying reversible particle model has, at 
equilibrium, phase transition for /3 > /3 C = see [22]. 
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We denote for / : Sn — »• R, x, y e Ayv and 77 e Sat, 

The microscopic dynamics is specified by a continuous time Markov chain on the 
state space Sn with infinitesimal generator given by 

Cn = L^jv + L-^jv + L+,jv , (2.3) 

where for function / : Sn — > R 

(L^/) fa) = 2 Cftfol/J^KV^/Jfo)] , (2.4) 

xeA N ,y£A N 

with rate of exchange occupancies C^- given by 

<( a; ,y;r 7 )=e X p|-^[ J ff J v(^)-J?iv(77)]} , (2.5) 
where iJjv is the Hamiltonian (2.2); 

(L-,iv/)(r?) = c-C^-JV))^ - ^) "/(»/)]. 
(L +>iV /)(r ? ) = C+ (r / (7V))[/(a A, 7 ? )-/(7 ? )], 

where for any p± e (0, 1), c± : {0, 1} — >■ R are given by 

c±(C) := P±(1-C) + (1-P±)C- 

The generator L^jv describes the bulk dynamics which preserves the total number 
of particles, the so called Kawasaki dynamics, whereas JL^jv, which is a generator 
of a birth and death process acting on Tjv, models the particle reservoir at the 
boundary of A^r. The rate of the bulk dynamics {C N (x, y; rf), x £ An, y £ An}, 
see (2.5), satisfies the detailed balance with respect to the Gibbs measure associated 
to (2.2) with chemical potential A € R 

^n X (v) = — ^ cxp{-f3H N (v) + A V{x)} , V G S N , (2.6) 

^jv xeAjv 

where Z N ' X is the normalization constant. For the bulk rates this means 
C N (x, y, V) = e-^ H ^- H ^C N (y, x; rf >*<). 

For the generator it means that L^g jy is self-adjoint w.r.t. the Gibbs measure (2.6), 
for any A € R. The corresponding process is denoted reversible. 

For a positive function p : An — > (0, 1) we denote by v N (-) the Bernoulli measure 
with marginals 

v$\r ) {x)=l)=p(x/N) = l-v N {r 1 {x) = Q) , x G A N , r, e S N . (2.7) 

The Bernoulli measure vp~ ( v p ~ ) is reversible with respect to the boundary gen- 
erator L +i tv (L-,jv)' 

The Markov process associated to the generator Cn, see (2.3), is irreducible and 
we denote by p N at = p N at ((3, p-,p + ) the unique invariant measure. In the notation 
we stress only the dependence on the parameters relevant to us. This means that 
for any / : Sn — > R 

r £ N f(v)dpT t (v)=0, 
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but the generator Ln is not self-adjoint with respect to p% at . The corresponding 
process is called non reversible. The only case where the process is reversible is 
when p_ = p + and (3 = 0. In such case the product measure associated to p_ = p + 
is invariant and the process with generator Cn is also reversible. 
We denote by M. the space of positive densities bounded by 1: 

M:= {p€ -Loo([-l, l],du) : 0< P <1} , 

where du stands for the integration with respect to the Lcbesgue measure on [—1,1]. 
We equip M. with the topology induced by the weak convergence of measures and 
denote by (•, •) the duality mapping. A sequence {p n } C M converges to p in M if 
and only if 

(p n ,G) = [ p n (u)G{u)du^ (p,G) 

J A 

for any continuous function G : [— 1, 1] —¥ K. Note that Ai is a compact Polish space 
that we consider endowed with the corresponding Borcl cr-algcbra. The empirical 
density of the configuration r/ G Sn is defined as ir N (77) where the map n N : Sn — > 
A4 is given by 

JV-l 

n N (r ] )(u) := £ r,(x) l{ [f - ^, § + ^) }(«) , 

in which 1{A} stands for the indicator function of the set A. Let {rj N } be a sequence 
of configurations with -q N G Sat. If the sequence {tt^^)} C M. converges to p in 
M. as iV — > 00, we say that {r/^} is associated to the macroscopic density profile 

2.3. Functional Spaces. Fix a positive time T. Let D([0,T],A^), respectively 
D([0, T] , 6>at), be the set of right continuous with left limits trajectories n : [0, T] — > 
A4, resp. (vt)t£[o,T] '■ [0,T] — > Sn, endowed with the Skorohod topology and 
equipped with its Borel a— algebra. Take pn on Sn and denote by (?7i)te[o,T] the 
Markov process with generator N 2 Cn starting, at time t = 0, by 770 distributed ac- 
cording to pn- Notice that the generator of the process has been speeded up by N 2 . 
This corresponds to the diffusive scaling. Denote by the probability measure 
on the path space -D([0, T], Sn) corresponding to the Markov process (%)te[o,Tl an d 
by the expectation with respect to P^^f • When pn = S n for some configura- 
tion 7j G S N , we write simply = P P S ' N and E^ N = E s f. We denote by tt n 

the map from D([0,T],S N ) to D([0,T],M) defined by n N (r].)t = n N {i lt ) and by 
QP£ = o (tt^)" 1 the law of the process (n N (rj t )) t>Q . For m G Loo([-1, 1]) 
and u G A we denote 

(J ncum *m)( M ) = / J ncum (u,v)m(v)dv. 
Ja 

We need some more notations. For integers n and m we denote by C n,m ([0,T] x 
[— 1, 1]) the space of functions G = Gt(u) : [0,T] x [—1, 1] — > R with n derivatives 
in time and m derivatives in space which are continuous up to the boundary. We 
denote by C*^ m ([0,T] x [-1,1]) the subset of C"' Tn ([0,T] x [-1,1]) of functions 
vanishing at the boundary of [-1, 1], i.e. G t (-1) = G t (l) = for t G [0,T]. We 
denote by C?- Tn ([0,T] x (-1,1)) the subset of C"< m ([0,T] x (-1,1)) of functions 
with compact support in [0,T] x (—1, 1). 
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2.4. Results. We denote x(p) = p(l _ p), a{p) = %x(p) an d W> res P- A/, tne 



gradient, resp. the laplacian, with respect to u of a function /. For G £ Gq' 2 ([0, T] 
[-1,1]), p 6 £>([0,T],.M) denote 

4(p, p ) := (pt, G T ) - (po, Go) - [ dt (pudtGt) 



x 



-/ dt(p t ,AGt) + p+ [ dtWG t (l)-p- f dtVGt(-l) 
Jo Jo Jo 

- 1 / V(pt), ( VG *) • v ( jncum * ft))*- 

Denote by A[ 0<T ] C L>([0,T];7W) the set 



(2.8) 



), T] = {peL 2 ([0,T],i/ 1 (A)) : VGeC 1,2 ([0,T]xA), 4(p,p )=o}. 



Theorem 2.1. For any sequence of initial probability measures (pn)n>1; the se- 
quence of probability measures (Q^^ v )at>i is weakly relatively compact and all its 
converging subsequences converge to some limit Q&>* that is concentrated on the ab- 
solutely continuous measures whose density p € A[o.t]- Moreover, if for any S > 
and for any continuous function G : [— 1, 1] — > R 

A^V^d^ £ n{x)G{x/N) - / 7 (u)G(u)du| > 5} = 0, 

/or aw initial profile 7 : A — > (0, 1), i/ien £/ie sequence of probability measures 
(Q^'n )n>i converges to the Dirac measure concentrated on the unique weak solution 
of the following boundary value problem on (t, u) £ (0, T) x A 

'd tPt (u)+ j 9V. { Pt (u)(l - Pt (u))V(J neum * Pt )(u)} = Ap t (u) 
< pt(Tl)= P T forO<t<T , (2.9) 
po{u) = 7(0) . 

Remark 2.2. ffy weofc solution of the boundary value problem (2.9) we mean 
4( 



4(ft l) = for G e C^dO, T] x [-1, 1]). 



Theorem 2.3. There exists f3o depending on A and J neum so that, for any (3 < /3q, 
/or any G G Cq([— 1, 1]), for any S > 0, 

lim p% at \\(ir N ,G) -(p,G)\ >S] =0, 

TV— >-oo L 1 

where p is the unique weak solution of the following boundary value problem 
'Ap(u)~[3V-{p(u)(l-p(u))V(J neum *p)(u)}=0, ueA, 
p(Tl) = Pt ■ 

We prove Theorem 2.1 and Theorem 2.3 in Section 4. Recall that the stationary 
measure p^°* depends on (3. 

Next we state the large deviation principle associated to the law of large numbers 
stated in Theorem 2.1. Let = j£ G : D([0,T],M) — > 1 be the functional 



8 



M. MOURRAGUI AND E. ORLANDI 



given by 

J§(n) := ^(Tr, 7) - I £dt (a(n t ), (VG t ) 2 ) , (2.11) 
and : D([Q,T],M) [0, +oo] the functional denned by 

/|(7r| 7 ) := sup J|(tt) . (2.12) 

GeC^ 2 ([0,T]x[-l,l]) 

To define the large deviation rate functional, we introduce the energy functional 
Q : D([0,T],M) -S- [0,+oo] given by 

Q(ir) = sup[ £ dt(n t ,VG t )-^£ dt(a(ir t )G t ,G t )} , (2.13) 

where the supremum is carried over all G £ C^°([0, T]x (— 1, 1)). From the concavity 
of cr(-) it follows immediately that Q is convex and therefore lower semicontinuous. 
Moreover Q(ir) is finite if and only if it e L 2 ([0, T];H 1 (A)) , and 

aw = 1 r« (2.i4) 

2J y_x CT(7T t (u)) 

If (2.14) holds, then an integration by parts and Schwarz inequality imply that 
(2.13) is finite. The converse needs to be proven, for a proof of it we refer to [4], 
Subsection 4.1. The rate functional I^i'h) '■ D([0,T],M) — > [0, oo] is given by 

-fH^m = s . (2-15) 

I +oo otherwise . 

We show in Lemma 5.6 that I^(tt\j) = if and only if 7r t (-) solves the problem 
(2.9) with initial datum ttq(-) = "/{■)■ 

We have the following dynamical large deviation principle. 

Theorem 2.4. Fix T > and an initial profile 7 in M. Consider a sequence 
{rj N : N > 1} of configurations associated to 7. Then, the sequence of probability 
measures {Q^'J? : N > 1} on D([0,T], A4) satisfies a large deviation principle with 
speed N and rate function defined in (2.15): 

151 llogQ^VeC) < -inf4H7) 

lim IbgQ^^eO) > -inf^H 7 ), 



N-too 



for any closed set C C -D([0, T], M.) and open set O C D([0, T], M). The functional 
^t('It) * s l° wer semi-continuos and has compact level sets. 

We prove Theorem 2.4 in Section 5. 



3. Basic estimate 



Next lemma states some properties of the potential J noum (-, •) easily obtained 
by its definition. 
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Lemma 3.1. The potential J nenm (. ; .) i s a symmetric probability kernel. Moreover 
for any regular function G : A — > R 7 we have the following: 

(y J neum (u,v)G(v)dv^j < J J neum (u,v)\VG(v)\dv. (3.1) 

Proof. The symmetry of J ncum follows immediately by the one of J. We have 

J noum (u, v) = J(0, v - u) + J(0, 2 - (u + «)) + J(0, 2 + (« + «)) , 

which is symmetric in u and u. We now prove that J neum is probability kernel. 
Fix u £ A, by a change of variables, 

/■(l-u)Al /■(3-«)Al 

J neun >,u)d« = / J(0,u)rfw+ / J(0,u)cfo 

/A J(-l-u)v(-l) 7(1— u)V(-l) 

r(-l— u)Al 

J(0, u)cfo . 

'(-3-u)V(-l) 

Suppose first that «£ [0,1], then 

y>l— ti />1 

J neun >,u)d« = / J(0,v)dv + / J(0,w)dw 



A J-l Jl-u 

1 

J{Q,v)dv = 1. 

-l 

and thus, J ncmn (it, •) is a probability. The proof for u € [—1, 0] is similar. It remains 
to prove (3.1): 

r J noum (u, v)G(v)dv) = J d u J ncum {u lV )G(v)dv 

J(u, v) - J(u, 2-v)- J{u, -2 - v) VG{v)dv . 
The result follows from the following incquatity, 

J(u, v) - [J(u, 2-v) + J(u, -2 - «)] < J neum (u, v) 
for all u, v e A. □ 

For any G:A4l and x, x + 1 € Ajv denote by V N G(- N -) the discrete gradient: 

X7 N G{x/N) = N[G((x + 1)/N) - G(x/N)] . (3.2) 

Next, we show that the rate C N of L^jv is a perturbation of the rate of the 
symmetric simple exclusion generator. 

Lemma 3.2. For any x G An , with x ± 1 € Ajv and 7/ € iS/v, 

c£(x, x± 1; tj) = 1 T ^Hx + 1) - nixflN- 1 ?* [( J neum ) **(»?)] (x/iV) + 0(iV- 2 ) . 

Proof. By definition of -ff/v, for all x, y € A^r and r] £ Sn, 

(V x , y H N )( V ) = ±{r)(x) - ry^)) 2 ^— (Jl, |) - J» cum (0,0)) 



+ (^)-r ? ( 2/ ))l ^ ^)[r™(^)-r™(|4)]. 

zeAjv 

This concludes the proof. □ 
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We start recalling the definitions of relative entropy and Dirichlet form, that are 
the main tools in the [11] approach. Let h : A — >■ (0, 1) and v£> ' be the product 
Bernoulli measure defined in (2.7). Given fi, a probability measure on Sn, denote 
by H(fj,\i/^ ') the relative entropy of fi with respect to ': 

ff(Ml4 (0 ) = su p{ J fivMdv)- log J e^pf\dr,)} , 

where the supremum is carried over all bounded functions on Sn- Since i/£ ^ gives 
a positive probability to each configuration, fi is absolutely continuous with respect 
to v N ^ and we have an explicit formula for the entropy: 

H(^) = [ logj-^W (3-3) 
J dv N y J 

Further, since there is at most one particle per site, there exists a constant C, that 
depends only on /i(-), such that 

HWT) < CN (3.4) 
for all probability measures fi on Sn (cf. comments following Remark V.5.6 in [15]). 

3.1. Dirichlet form estimates. One of the main step for deriving the hydro- 
dynamic limit and the large deviations, is a super exponential estimate which al- 
lows the replacement of local functions by functionals of the empirical density. 
One needs to estimate expression such as (Z, in terms of Dirichlet form 

N 2 (—£Nx/f{i]), \J f{v))p, N > where Z is a local function and (-,-) m jv represents 
the inner product with respect to some probability measure p N . In the context 
of boundary driven process, the fact that the invariant measure is not explicitly 
known introduces a technical difficulty. We fix as reference measure a product 
measure . see (2.7), where 9 is a smooth function with the only requirement 

that #(t1) = Pt- There is therefore no reasons for N 2 (—Cn\/J(j]), \J f {v)) v e (■) 
to be positive. Next lemma estimates this quantity. 

Define the following functionals from h € L 2 (v) to R + : 

N-l 



V , N (K v)=Y,j (HV X ' X+1 ) - Hv)) 2 dv{rf) , 

x=-JV' 

V + . N {h, v) = \j c + ( V (N)) (h^-^) h{ V )) 2 du( V ) , ^ 3 - 5 ) 



P-,w(M = 2 / C -(V(-N) {h{a- N+1 r 1 )-h{r 1 )) 2 dv{r ] ). 

Lemma 3.3. Let 6 : A —> (0, 1) be a smooth function such that 0(^pl) = p T . 
There exists a positive constant Cq = Co(||V#||oo) so that for any a > and for 

[ m^mduf\r,) < -(1 - a)Vo,N(f,^ } ) + —N- l \\ff H . ) t (3.6) 
Js N a { N 1 
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Proof. The proof of (3.7) follows from the reversibility of the Bernoulli measure 



/(•) 



' with respect to L± jy. Next, we show (3.6). By Lemma 3.2, 



N-l 



=-JV 



N 



N-l 

E 

x=-N 



(3.8) 



JX,X-\-l 



/)(»?) mdvViji) 



for some positive constant Ai depending only on /3 and J. We write the first term 
of the right hand side of (3.8) as 

N-l 



E 

x=-N 
N-l 



jx,x+l 



/)(») d v %\rj) 



(3.9) 



x=-N 



E RN(x,x + l;9,r,) (V*'* +1 /)M f{^ +1 )du^'\ri) 



where 



ie JV (x,x+l;e,T?)= [i-e-^lW)^))] , (3.10) 
A is the chemical potential defined by 

A(r)=log[r/(l-r)] (3.11) 
and S7 N stands for the discrete derivative defined in (3.2). By the inequality 



a -2 1 B 2 



for all A, B € K and a >0 , AB < -A 

2 2a 

and Taylor expansion, the formula (3.9) is bounded by 

N-i „ „ . 

A 2 



(3.12) 



l-~) E / d^\ V ) + ^N-^f\\l 2{vr) (3.13) 



for all a > 0. Here A 2 is a positive constant. 

The second term on the right hand side of (3.8) is handled in the identical way. 
It is bounded by 

N-l 



-N 



E / \(V x ' x+1 mYdv%\v) + ^N-i\\f\\< [ H 



The lemma follows from (3.8), (3.9), (3.13) and (3.14). 
Denote for h € L 2 (u) 



(3.14) 

□ 



iV-l 



Vf,, N (h,u)= E / C p N (x, x + 1; i]) (h{if> x+1 ) - h(ji)) dv(ri) . 



x = -N 

Lemma 3.4. There exists a positive constant C\ = C±(0,J) such that, for any 
measure v and for h £ L (y), 

(1--^-)2?o,jv(M) <T>p, N (h,u) < (l + -±)V 0ttf (h,v). 

Proof. The proof is elementary since |C^-(x, x + 1, rf) — l| is uniformly bounded in 
N, x and r\. □ 
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Lemma 3.5. Let p, po : A — > (0, 1) be two smooth functions. There exists a positive 
constant C = Cg(|| Vpolloo, ||V/?||oo) such that for any probability measure p N on 



dp 



N 



dp 



/>(•) 



/>(■) 
77 > N 



< 2 2), 



(3.15) 



P0(") 



('/I 



< 2 



TV 

Proof. Denote by /fa) = j^{ V ) and fcfa) = -fefa). Since f( V ) = h(r,)% 

au N ai/ N 

we obtain for — N < x < N — 1 the following 

\/M ? 7^ +1 ).R 2 (:e,:e + 1; 77) + Vx,x+iVM»?)] % dv p ^\rf) 
V x , x+1 Vh(r,)} 2 dis p N o{ -\n) 



+ 2 / fc(»f' a!+1 )[i2 JV (x > x + l;p, )7 )] 2 di/« ( - ) (»7) 1 

where 

J2a(x,i + 1; 77) = exp {(l/2)iV- 1 V Ar [A(p(x/7V)) - \{p (x/N))]V x , x+ iv(x)} - 1 

and A is the chemical potential defined by (3.11). We conclude the proof using 
Taylor expansion and integration by parts. □ 

3.2. Superexponential estimates. For a positive integer I and x £ An denote 

A e {x) = A NJ {x) = {y£A N : \y - x\ < £} . 
When x = 0, we shall denote A^(0) simply by A^, that is, for all 1 < t < N, 

At = Aiv^(0) = {—I, ■ ■ ■ ,£}. 
Denote the empirical mean density on the box An(x) by rj l (x): 



1 



IM*)I 



(3.16) 



y£A e (x) 



For a cylinder function 'J, that is a function on {0, 1} Z depending on r)(x), x £ Z, 
only trough finitely many x, denote by ^(p) the expectation of ^ with respect to 
v p , the Bernoulli product measure with density p: 



*(p) = £<[*] . 
Further, denote for G £ C([0,T] x [-1, 1]) and e > 

1 



(3.17) 



(3.18) 



where the sum is carried over all x such that the support of t x \P belongs to An and 
[•] denotes the lower integer part. 



KAWASAKI PROCESS WITH NEUMANN KAC INTERACTION 



13 



Proposition 3.6. Let {jUjv ■ N > 1} be a sequence of probability measures on Sn- 
For every 5 > 0, 



1 



lim lim — logP^lf 

_ , n at AT ° MiV 



> 5 



-00 . 



Proof. Fix c > that will decreases to after e and a smooth function p c : A — > 
(0,1) which is constant in An_ c ) = [— 1 + c, 1 — c] and equal to p± at the boundary, 
i.e p c (±l) = p±. The constant can be arbitrarily chosen and we denote it 70. 
Divide Ajy in two subsets, Ar(i_2 C )jv] an d A^ \ h-u\-2 C )N] an d split the sum over x 
in the definition of V N ' into the sum over these two sets. Since 



sup 

V7,E,JV,xgAj 



we have that 



E 

3.€A N \A[ (1 _2c)iV] 



{g s (x/N) [t x V( V ) - *(r,[ £jv ](x))] } 



dsG s (x/N) [ Tx *(f/) - *(?7 [eArl (x)) 



< 00 , 



< cTK (3.19) 



for some positive constant Kq which depends on G. By Chebyshev exponential 
inequality, for all a > 0, 



1 

N 



V N.f( s ' 1 ls)ds 



>6 



<-a(5- TK c) 



N 



log 



t/3, JV 



exp aAT 



where 

= ^ E g s {x/n) [ Tx *{n) - *(v sN (x)) 

i£A[(i_ 2c )jv] 

It is immediate to see that the Radon-Nikodym derivative 



(3.20) 
(3.21) 



(7;/ 



Pc(-) 

TV 



< gA^i(c) 



for some positive -f^i(c) that depends on c. The right hand side of (3.20) is bounded 

by 



1 



-a(S-TK c) + Ktic) + — logP^; 



exp aN 



(3.22) 



Since e' x ' < e x + e x and 

HmA''" 1 log{ajv + bjv} < max{hmiV _1 loga^v , En A^ 1 log b N }, (3.23) 

we may remove the absolute value in the third term of (3.22), provided our estimates 
remain in force if we replace G by — G. Denote by 

(£n) s = 7j {£n + £*n) 
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where C* N is the adjoint of £jv in L 2 (v p ^ ^ By the Feynman-Kac formula, 



— logP , 

N 8 < ( - 



T 



1 



expjaiV / V^ , "(«,7? a )d8j < -^J X N , e (G s )ds, (3.24) 



'i' 



TV 

( 

In this formula the supremum is carried over all densities / with respect to v p ^ ' 

2gative contribution, it is enough, to g( 
result, to choose c such that c < jpfr- and to show that, there exists M > that 



where \ N ^(G S ) is the largest eigenvalue of {N 2 (C N ) S + Na¥ c ^ 9 (G s , •)}• By the 
variational formula for the largest eigenvalue, for each s € [0,T], 

^Aiv >£ (G s )=sup{ f aW%^{G,,T,)f(T,y/-\dn) + N(C N ^f, ^f)^,} • 

(3.25) 

V N ■ 

By Lemma 3.3, since (3.7) gives a negative contribution, it is enough, to get the 
result, to choose c such that c < and to sho 
depends only on G and c, such that, for all a > 

m En supf / a V^*(G s ,7 ? )/( ?? X c() (^) - iVA^CV/,^ )) < A/ - 

We then let a | oo. Notice that for N large enough the function Vjy £ ' (G s ,rf) 
depends on the configuration r\ only through the variables {r](x), x e A(i_ c )at}. 
Since p c is equal to 70, in A c , we replace v p ^ ' in the previous formula by v]^ 
with respect to which the operator I^o,n is reversible. Therefore 2?o,iv( - > v~n) IS the 
Dirichlet form associated to the generator Lo^at. Since the Dirichlet form is convex, 
it remains to show that 

Bm HE supf f aY&"(G a ,ri)f( V )v%(drj) - NV , N (Vf, "£ (0 )} = 0, 

e— »0 A/— ¥00 f I J ' 

for any a > 0. This follows from the usual one block and two blocks estimates (cf. 
Chap 5 of [15]). □ 

For x = ±iV, a configuration r\ and £ > 1, let 

W^fa) = |^(±JV)-p±|, (3.26) 

where, see (3.16), rf(N) = -^f{v(N -£) + ■■■ +r)(N)} and rf(-N) = I ^ I {r)(-N + 
£) + ... +V (-N)}. 

Proposition 3.7. Fix a sequence {/in ■ N > 1} 0} probability measures on Sm- 
For every 5 > 0, 

-00 . 



lim lim — logPf' W 







Proof. Consider first the limit with the term W^' . Fix a smooth function 7 : 
A — >■ (0, 1) such that 7(— 1) = and 7(1*) = p + for u S [0, 1]. Since the Radon- 
Nikodym derivative % is bounded by exp(iV-K'i) for some positive constant K\, 

dv N 

it is enough to show that 
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We follow the same steps as in Proposition 3.6. Applying Chebyshev exponential 
inequality and Feynman-Kac formula, the expression in the last limit is bounded 
for all a > by 

-aS + ^ V e (a) , (3.27) 
where for all a > 0, i Ajv !£ (a) is the largest eigenvalue of the vjj- -reversible operator 



f^N(C N ) s (f) + a(w+ J )f. 



Here (Cn) s is the symmetric part of the operator Cm in L 2 (vjj-^). By the variational 
formula for the largest eigenvalue, we have 



l~\ N Aa) = sup! / aWp'i^fivy^idrj) 



N 



N ) 



In this formula the suprcmum is carried over all densities / with respect to v]j ■* . 
By Lemma 3.3, we just need to show that, there exists M > 0, such that, for all 
a > 

< M 



lim lim sup<^ / aW N ' {v)f(v)^N ( dr l) 



~ NV 0tN (y/f, vl ( - ] )-NV +tN (^]y N { - ] )) 

Recall that the profile 7 is constant equal to p + on [0, 1]. Since W^' (77) depends 

y]j ■* in the previous formula by v 9 ^ 



only on coordinates in a box Af(N), we replace vjj ■* in the previous formula by u 9+ 



On the other hand, v 9 ^ is reversible for Lo.jv + ^+,n and therefore 2?o,jv( - , v^n) + 
T> +: n(- , v p ^) is the Dirichlet form associated to the generator L 0i w + L+,jv- Since 
the Dirichlet form is convex, it remains to show that 

1ST EE su P { f aWp £ ( v )f( v y N +(d v ) 

I— >oo iV— ^oo j I J 

- NV , N {^f,v p +)-NV + , N {^]y N + )} = 0. 

for any a > 0. This follows from the law of large numbers by applying the same 
device used in the proof of the one block and two blocks estimates, (cf. Chap 5 of 
[15], and Lemma 3.12, Lemma 3.13 in [20]). □ 

3.3. Energy estimate. We prove in this subsection an energy estimate which is 
one of the main ingredient in the proof of large deviations and hydrodynamic limit. 
It allows to prove Lemma 4.3 and to exclude paths with infinite energy in the large 
deviation regime. For 8 > 0, G G C£°([0,T] x A) define 

SgW = / dt{n t ,VG t ) - 5 dt{a(n t )G t ,G t ), (3.28) 
Jo Jo 

Q 5 (n) = sup few} ■ (3-29) 

G£C»([l),T]xA) L J 



Notice that 



where Q(-) is defined in (2.13). 



Q s (n) = lew, 
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For a function m in Ai, let m £ : A 



be given by 



™ £ 0) = i I 

zt J\u—e,u 



m(v) dv 



+e]HA 

When u G [— 1 + e, 1 — e], rn e {u) = (to * i e )(it), where t e is the approximation of the 
identity defined by 

L e{u) = ^ !{[-£, e]}(u) . 

Lemma 3.8. There exists a positive constant C± depending only on p± so that for 
any given Sq > 0, for any S, < 5 < 5o, for any sequence {i] N G Sjy : AT > 1} and 
for any G G C r °°([0,T] x A), we have 



lim lim 



1 



e->0 7V^oo N 



loaP 



fi,N 



exp (8 NQq(tt n * t e )) < Gi(T + l) 



Proof. Assume without loss of generality that e is small enough so that the support 
of G(-, •) is contained in [0, T] x [—1 + e, 1 — e]. Let 6* : A — > (0, 1) be a smooth 
function such that 6(=pl) = p T . Since v^j'\ri N ) > exp{— C(Af} for some finite 
constant C[ depending only on 9, it is enough to prove the lemma with 



place of 

Set #1(77) = [r/(l) - r;(0)] 2 and note that $ x (o) = = cr(a) = 2a(l - a), 

where v a is the Bernoulli measure with parameter a G [0, 1]. Denote B N ' e g the 
set of trajectories (?7i)tg[o,T] so that 

r-T 



G,*i 



[r?. G D([0,T],«Siv) : | ^ V^fo 7fc)dt| < ^} , 



where V^'* 1 is defined in (3.18). By (3.23) and Proposition 3.6, it is enough to 
show 



lim lim — log P' 3 ^ 



HB^} g } < d(T+l) 



Recalling the definition Q 6 G , see (3.28), we have 



rp rp N — l 

[ dt(7rf *t e ,VG) = / dt V {^(^-^(x + ^IGt^/^ + Oc^). 
Jo Jo __ 



a; = -Ar+l 



Further on the set -B 



G.*i 



Jo / di(ff(7rf * O 



7 1 AT— 1 

,G 2 t )>5 Q f dt- G^x/N)^^) 

J „ AT" | 1 



x=-iV+l 



-<5 O G2 (jV j£ )--, 

where Oq(e) is absolutely bounded by a constant which vanishes as e | and 
Oq2(N, e) is is absolutely bounded by a constant which vanishes as A~ f oo. There- 
fore to conclude the proof it is enough to show that 



lim — Ioe 

N^oo N 



$(■) 



cxp^N j dtV£(t,r) t ) 



< CiT 



(3.30) 
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for any 5 < So, where 

N-l r2 W-i 

V*(t, V ) = 6 J2 G t (x/N)[ v (x)-r,(x + l)}-^ J2 GHx/N)^^). 

x = -N+l x=—N+l 

Now, observe that Vq = Vg G . Therefore, to prove the lemma, we need to show that 
for any smooth function G, 

EE ilogP^fexpfiV / dtV^rjt))} < C X T (3.31) 

for some constant C\ that not depends on G. By Feynman-Kac formula and the 
same arguments used in the proof of Proposition 3.6, the expression of the limit in 
the right hand side of (3.31) is bounded above by 

cftsupf / V&(t t ri)f a (<n)u#\dri) + N(C N fJ) i/( ,\ , 
o / W N > 

where the supremum is over all functions / in L 2 (i/^ ') such that (/, f) v e{) = 1. 

By Lemma 3.3, we replace N{CnI, f) H ■) by — N(l — 6)2?o.A'(/, v k ) + tK where 
b € (0,1) is arbitrarily chosen and Co is a constant depending only on p±. It 
remains, therefore, to show that 



lim / dtsupf / V£(t,r))f(r,yV(dr,) - N(l - b)V . N (f, v 9 ") } < C X T, 

(3.32) 



We split 



where 



V£(t,r,)f 2 ( V )v 9 N { - ) (d V )=I 1 -I 2 , 



N-l . 



x=-N+l 

1 N_1 r 

x=-N+l J 

We estimate Ii in term of J2 and 2?o ,jv(/> ^n'^)- By changing variables 7/ = 7y a: ' a:: + 1 ; 
we have that 

AT-l . 

^i=2 E G t (x/iV) /{^(.t)-^ + 1)}{/ 2 W-/ 2 (^ +1 )}< ( - ) W 

E=-JV+1 



+ | E G *(*/A0 fMx)-v(x + l)}RN(x,x + l ] e,7 ] )f(7 ] )dv e N { -\r l ), 

(3.33) 



a;=-AT+l 
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where Rn(x, x+1; 9, -q) is defined in (3.10). For the first term of (3.33), by inequality 
(3.12) and Taylor expansion, we have 



JV-l 



\ J2 G t {x/N) f{ V {x)- n {x + l)}{f\r 1 )-f\r ] ^ +1 )}dvf\r 1 ) 



x=-N+l 

< aN T) (fJth 

N-1 



l — G t( x / N ) [r^i(v)[m + f(v x ' x+1 )] 2 d^\v) (3.34) 

AT" I 1 *> 



4aN 

x=-N+l 

<^ , w (/,^)) + -Lc (G ) 

N-l 



± J2 G 2 t (x/N) J 'rMv)f 2 (v)du e N { -\v) 



x=—N+l 



where C(G) is some constant that depends on G. For the second term of (3.33), 
by (3.12) and Taylor expanding Rn we have that 



N-2 



x=-N+l 



\ G t (x/N) f{r ] (x)-r,(x + l)}R N (x,x + l- > e,r ] )f(r,)d^\r 1 ) 

+i J 

i N ~ 1 r 

— Yl J G t {x/NfT x Mri)f{v)dv e <;-\ri), 



<Ca + 



x=-N+l ' 



1 



Cai — 12 
a 



(3.35) 

for all a > 0, for some positive constant C depending only on p±. Taking into 
account (3.33), (3.35) and (3.34) we have 



2 r aAT , , 
h < -h + --V . N (f 
a 4 



We conclude the proof, by taking a — 2 and b = i in (3.32). 
The following corollary allows to show Lemma 4.3. 



(3.36) 
□ 



Corollary 3.9. Fix a sequence {Gj : j > 1} C C£°([0,T] x A), S Q > and a 

sequence {n N £ Sn ■ N > 1} of configurations. There exists a positive constant C\ 
depending only on the values p^, such that for any < 5 < 5q and any k > 1 



lim lim — log 



cxp(<5iV sup Qg° (tt n * i e )) < C x (T + l). (3.37) 

i<j<k 3 J 



e-rt) JV->oo TV ^ 

Proof. From (3.23), the limit in (3.37) is bounded above by 

max ( lim lim" — logP^f [exp (S NQ S J1 (n N * u))] ) 



By Lemma 3.8 the thesis follows. 



□ 
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4. HYDRODYNAMIC AND HYDROSTATIC LIMITS 

Wc prove in this section the hydrodynamic and hydrostatic limit for our system. 
The proof is based on the method introduced in [11] for the hydrodynamic limit 
and in [7] for hydrostatic, taking into account, as explained in the introduction, the 
lack of comparison and maximum principle of (2.9). 

4.1. The steps to prove Theorem 2.1. Following [11] we divide the proof of 
the hydrodynamic behavior in three steps: tightness of the measures (Q^), an 
energy estimate to provide the needed regularity for functions in the support of any 
limit point of the sequence (Qu'^)> an( i identification of the support of limit point 
of the sequence (Q^) as weak solution of (2.9). We then refer to [16], Chapter 4, 
that present arguments, by now standard, to deduce the hydrodynamic behavior of 
the empirical measures from the preceding results and the uniqueness of the weak 
solution to equation (2.9). 

Lemma 4.1. (Tightness) The sequence {Q^) is tight and all its limit points 
Q^'* are concentrated on the following set 

Q"'*{tt : < 7r t («) < 1, te[0,T], ue[-l,l]} = l. (4.1) 

Wc then show that Q 13 '* is supported on densities p that satisfy (2.9) in the weak 
sense. 

We start defining for G G C 1,2 ([0,T] x A) and e > 
Bf' N = [ G T (u)7r N (r lT )(u)du- f G (u)n N (r l0 )(u)du 

J A J A 

r T 



d s G s (u)TT (r] s )(u)duds — / / AG s (u)tt (n s )(u)duds 

Jo J A 

E / T (VG s )(.x/iV){ ( T(^](.))V w (J— *n N (r, s ))(x/N)}ds 

x<£A N J ° 

T 

dt[p + (VG t )(+l)-p-(VG t )(-l)] , 

o 

(4-2) 

where r)^ N \x) is the local mean defined in (3.16) and \7 N G s (x/N) stands for the 
discrete gradient of G s (x/N) defined in (3.2). 

Proposition 4.2. (Identification of the limit equation). For any function G 
in Co '■ ([0,T] x A) and any 8 > we have 

En IW F?,' N (\B?' N \ >5) =0. (4.3) 

The last statement is an energy estimate. Every limit point Q^'* of the sequence 
(Q^) is concentrated on paths whose densities p € £ 2 (0, T; iJ 1 (A)). 

Lemma 4.3. (Energy estimate) Let Q 13 '* be a limit point of the sequence (Q^)- 
Then, 



Q 



13,* 



L 2 (0,T;H\A)) =1. (4.4) 
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4.2. Proof of Proposition 4.2. Let * be a limit point of the sequence (Q }f ) 

and assume, without loss of generality, that Q^'n converges to Q@'*. Fix a function 

G in Cq' 2 ([0,T] x A). Consider the ¥^' N martingales with respect to the natural 

filtration associated with (r/ t )te[o,T], Aff = M^' N ' and TVf = A/f' N ' p , £ € [0,T], 
defined by 

M G =< Trf , G t > - < Trf , G > - / ( < Trf, 9 s G s > +iV 2 £^ < Trf , G s > )ds , 

Jo 

M G = (Mf) 2 

y {iV 2 £&( < Trf , G s > ) 2 - 2 < Trf , G s > 2V 2 £^ < ^ Gs >} ds . 

(4.5) 

A computation of the integral term of A/]f shows that the expectation of the qua- 
dratic variation of M t G vanishes as N f 0. Therefore, by Doob's inequality, for 
every S > 0, 



lim Pf' w 

A^OO 



sup |il/ t | > S = . (4.6) 

0<*<T 



Since for any s £ [0, T] the function G s vanishes at the boundary of A, a summation 
by parts permits to rewrite the integral term of M t as 

< Trf ,<9 S G S > ds 
ft N-i 

~ N { 12 (^ N Gs)(x/N)C^(x,x + l,r] s )(V x , x+1 ri s (x))jds, 

J ° x=-N+l 

where V N is defined in (3.2). 

From Lemma 3.2, a summation by parts and Taylor expansion permit to rewrite 
the last expression as 

ft rt 



0(N) + I < ?rf ,<9 S G S >ds+ f < Trf , AG S > ds 
Jo Jo 

+ f {-VG s (1)tj s (A0 + VG s (-l) Vs (-N)]ds 



+ 4v f { S (VG s )(x/iV)(V :c , :!;+1??s (x)) 2 V JV (J noum *Tr JV (r ?s ))( a ;/7V)}d S . 

Next, we use the replacement lemma stated in Proposition 3.6 and Proposition 3.7. 
We obtain that the integral term of the martingal M t G can be replaced by 



t 

< Trf , d a G a >ds+ I < Trf , AG S > ds 
Jo 

■J {-VG s (l)p + + VG a (-l)p_}ds 

•f /a7{ £ (VG s )(x/iV)a(^(x))V JV (J ncura *Tr w (r, s ))(.T/iV)} 



This concludes the proof of the lemma. □ 
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4.3. Steps to prove Theorem 2.3. Let p s ^ at = /i^ at (/3, p-, p+) be the unique 
stationary measure of the irreducible Markov process (r]t)t>o with generator Cn- 
From Tchebyshev's unequality, we need to show that 



lim E^ N 

N->oo 



(n N ,G)-{p,G) 



= 0, 



(4.7) 



where E^ N stands for the expectation with respect to the stationary measure 
H N at . It is enough to prove that any subsequence of the sequence of real numbers 
in the limit (4.7) vanishes. Without loss of generality we consider a sequence in 
(4.7) as a subsequence along which the limit exists. 

Denote by QP,^,stat . = qP.n^ the pro bability measure on the Skorohod space 

D([0,T];M) induced by the Markov process (tt^) = (ttjv (vt)) j when the initial 
measure is p% at . 

By the first part of Theorem 2.1 we have that all limit points of the sequence 
q/3, N.stat are concentrated on -A[o,t] f° r an y T > 0, i.e all its limit points are 
concentrated on the weak solutions of the hydrodynamic equation for some unknown 
initial profile. 

-sequence converging to a limit point which we denote 
by QP'*' stat . Note that different subsequences might have different limit points. 
Let j3 small enough and denote by p the unique stationary solution of (2.9), see 
Theorem 6.2. By stationarity we have for any S > 0, 



stat 

E tJ ' N k 



(n N ,G)~(p,G) 



" stat 



Since the integrand is bounded we have the following: 

-<P>G)\} 
E^ stat {(( PT7 G)-(p,G) 



lim E'\ 

k-toc Vn,. 



< ||G|| 2 ^" etet {| 
< ||G|| 2 e^) T 



Pt 



1 



{-4[0,T]} 



(p)} 



by Theorem 6.2 and ||u||2 denotes the L 2 (A) norm of v. Then letting T 
show the thesis. 



(4.8) 



oo we 
□ 



5. Large deviations 

In this section we prove some properties of the rate function and we present the 
main steps to derive the large deviations results. 

Let L 2 (A) be the Hilbert space of functions G : A — > R such that J A \G{u)\ 2 du < 
oo equipped with the inner product 



(G, J) 2 = / G(u) J(u)du . 
J A 

The norm of L 2 (A) is denoted by || • I2. 

Let -ff 1 (A) be the Sobolev space of functions G with generalized derivatives VG 
in L 2 (A). i? 1 (A) endowed with the scalar product (•, -)hi, defined by 

(G,J) H i = (G,J) 2 + (VG, VJ) 2 , 
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is a Hilbert space. The corresponding norm is denoted by || • ||#i . Denote by i?o(A) 
the closure of (A) (the set of infinitely diffcrentiable functions from A to K with 
compact support in A) in H 1 (A). Denote by H~ 1 (A) the Hilbert space, dual of 
i?o(A), equipped with the norm || • ||_i 



sup \2(v,G)- ltl - [ \VG(u)\ 2 du\ 

EC" (A) I J A J 



GeC»(A) 

where (v,G)-i t \ stands for the duality between Hq and H~ x . Fix T > 0. For a 
Banach space (B, || • ||b) we denote by L 2 ([0, T],B) the Banach space of measurable 
functions U : [0, T] -> B for which 

ll^lli^lCTLB) = / HCtllid* < 00 
Jo 

holds. 

5.1. Properties of the rate function. Denote 

={ir£D([0,T],M):ir (-) =l(-)l Jr t (±l) =p±,t G (O.T]}. 

Lemma 5.1. Let ir be a trajectory in D([0,T], A4) such that Jy(7r|7) < oo. T/ien 

Pi 
7 



7r belongs to B p± 



The proof is similar to the one of Lemma 3.5 in [2]. To prove the lower- 
semicontinuity of the rate function, we need the next results 

Lemma 5.2. For any (3 > ; there exists a constant Cq = Cq(/3) such that 
IM^ < C {^(7r|7) + fi(7r)}, Q(tt) < C {l + ^(tt| 7 )} 
for all 7T m D([0,T],M). 

Proof. The proof is the same as in Proposition 4.3. [21] or Theorem 3.3. in [19], or 
Lemma 4.9. in [4] □ 

Lemma 5.3. Let {p n : n > 1} be a sequence of functions in L 2 ([0,T] x A) such 
that 

dt\\pl l \\ 2 H i + I dtWd^W < c 





for some finite constant Cq and all n > 1. Suppose that the sequence p n converges 
weakly in L 2 ([0,T] x [—1,1]) to some p. Then, p n converges strongly in L 2 ([0,T] X 
[-1,1]) to p. 

Proof. Recall that H^A) C L 2 (A) C ^ _1 ( A )- % [ 25 > Theorem 21.A], the em- 
bedding -ff 1 (A) C L 2 (A) is compact. Hence, by [23, Lemma 4, Theorem 5], the 
sequence {p n : n > 1} is relatively compact in L 2 (0,T; L 2 (A)). In particular, weak 
convergence of the sequence {p n : n > 1} implies strong convergence. □ 



Theorem 5.4. The functional 7^(-|7) is lower semicontinuous and has compact 
level sets. 

Proof. Theorem 5.4 is proven applying Lemma 5.2 and Lemma 5.3. See Theorem 
3.4. in [19] or Lemma 4.2. in [4]. □ 
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We provide an explicit representation for the rate function iy('l'y) when it is finite. 



Fottt £ D([0,T],M), denote by H the Hilbert space induced by Cq' ([0,T] 
[— f , 1]) endowed with the inner product (■, Oo-fir) defined by 

(F,G) aM = [ dt(<j(n t ),X7F t -VG t ). 

Induced means that we first declare two functions F, G in Cq' 2 ([0, T] x [— 1, 1]) to 
be equivalent if (F — G,F — G)^) = an d then we complete the quotient space 
with respect to the inner product (•, •)„■(„•). The norm of iJg(cr(7r)) is denoted by 

II ' llcr(Tr)- 

Lemma 5.5. Take it £ D([0,T], M.) with I^(ir\j) < oo. Then, it is uniquely de- 
termined a function F in Hq(o~{tt)) such that ir is the weak solution of the following 
boundary value problem: 

Att- V- {cr(7r)[fV(J" etim *7r) + VF]} in Ax(0,T), 
7(0 <n A, (5.1) 
P± for 0<t<T. 




4(*h) = \\\F\\l M = \ £dt (a(n t )VF t ■ VF t ] 



(5.2) 



Proof. By assumption I^(tt\j) = i^(ir\j) < oo, defined in (2.12). Then one pro- 
ceeds as in [16] with the only difference that because the boundary conditions the 
space is Hq(o-(tt)). □ 



Lemma 5.6. Let p £ L 2 ([0,T], H 1 (A)) be the weak solution of the boundary value 
problem (2.9) then 

4ip\l) = It(p\7) = 0, and Q(p) < oo. 

Further if Ij,(p\j) = 0, then p £ L 2 ([Q,T], H 1 (A)) is the weak solution of the 
boundary value problem (2.9). 

Proof. We start showing that if p £ L ([0, T],iJ (A)) is the weak solution of the 
boundary value problem (2.9) then Q(p) < oo. Take F(p) = plogp+ (I — p) log(l — 
p), for p £ [0, 1]. Since J A F(p t (u))du is a bounded quantity for all t £ we have 
that 

I dt Jt f A F ^ u ^ du = - F(p Q (u))]du. 

Notice that 

F'(p)=log — ^ and F"{p)^—^- = ^- 
(1 - P) P(l - P) XKP) 

are not uniformly bounded for p £ (0, 1). Therefore we need some care to derive 
F(p t (u)) with respect to t. We consider a sequence of smooth functions 

F n (p) = (i + -y\p + -)iog(p + -) + (i-p+-)\og(i-p + -) 
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so that lim n _>. 00 F n (a) = F(a). We have 

T dtj t Jj n { Pt {u))du = dt Jj' n {p t {u))j t p t {u). (5.3) 

To avoid boundary terms, take a smooth function &(•) defined on a neighborhood 
of [-1, 1] such that 6(=pl) = p^ and < p~ < &(•) < p + < 1. Denote 

U n (t,u) = FMu)) - F^(b(u)). 

We have 

T dt I F' n { Pt {u))j t p t {u) 



J A 

T 



dt I Un {p t {u))-p t {u) + J dt J^(b(u))j t p t (u) (5.4) 

dt J U n (pt(u))—p t (u)du + J F^{b(u))[p T (u) - p Q (u)]du. 
Taking into account (5.3), (5.4) and U„ G L 2 ([0, T], J/q), we get 

T 



dtj t F n {p t {u))du - F' n {b{u)) [ PT {u) - P) («)] 

T dt / W n (i,«) [Vpt(u) - 0/>t(u)(l - p t (u))(J ncum * V Pt )(u)] 

J A 

Denote Xn{a) = (a + — )(1 + r — a)- We have that 

V/9 t (w) V6(«) 



VC/„(t,u) 



(5.5) 



Xn(Pt( w )) (&(«)) 

Taking this into account and collecting the above estimates, we obtain 

J A Xn{pt(u)) 

<- [ [F(p T (u))-F( Po (u))}du+ f F^b(u))[ PT (u) - Po (u)]du 

J A J A 

T , f Vb(u)-Vp t {u) , , x 

7a Xn(oW) 

+ /3 / dt ( Vp t (u)) * ( f ( "\\ (J ncum * V Pt )(")rf" 



- /3 / / V&H x( ^ ( "Pj J noum * V ft )(u)d« . 

JO JA Xn{b{u)) 

Since &(•) is bounded below by a strictly positive constant and above by a constant 
strictly smaller than 1, and since 

dt [ du\7p t {u){J ncum *\7p t ){u) <C 

J A 

for some constant C, we obtain, uniformly in n 

T dtf^^du<C 

J A XnKPtW) 
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for some finite constant C which depends only on b and T. To conclude the proof 
it remains to apply Fatou's Lemma and recall the definition of Q(p) given in (2.14). 
We have shown that Q(p) < oo. By Lemma 5.5 we conclude that (p|t) = 0. 
Similar arguments allow to prove the second statement of the lemma. □ 

5.2. Comparison between /£(-|7) and Next, we compare the rate func- 

tional Ij'(-\ r y) with the rate functional I^i'll) of the symmetric simple exclusion 
process (i.e. /? = 0). 

Lemma 5.7. For tt G D([0,T],A4), with finite energy Q(tt) < oo, we have 

2 ,r - ( 57 ) 



2 JV 16 7 



< 2/°(7r| 7 ) + dtj (V7r t ) 



Proof. Fix tt g D([0,T],M) with finite energy and G g C 1,2 ([0,T] x [-1,1]). 
Recall from (2.11), (2.8) and (2.12) the definitions of j£(ir), i% and I$(ir). By the 



inequality ab < ^a 2 + ^b 2 we obtain 



2 1 2 



4(^7) -^K7) 



/ di / a(^)(VG t )-V(J ncum *7r t ; 

JO J A 



<] f 1 dt f a(n t )(VG t ) 2 + Tdt [ a(7r t )[V(J ncum *7r t )] 2 . 

4 JO JA 4 JO J A 

Since for each u, J ncum (u,v)dv is a probability density on A and cr(-) < 1/2, by 
Lemma 3.1, Jensen inequality and Fubini's Theorem, 

?/ T ^j( ff(7rt)[V(jnCUm * 7rt)]2 

[ T dt I J neum *(Vn) 2 = ^- [ T dt [ (V7r t ) 2 . 

"JO J A 8 JO JA 



Hence 



< l° G W - \£ ' dt{a(n t ),(VG t ) 2 ) + ^ j\t jf (Vtt*) 2 
<i sup {2£°(tt)- J / dt(<r(7r t ),(VG t ) a )} 

Z GeCo ,2 ([0,T]xA) *>0 J 

= 2^(7T) + 4" / di / ( v ^) 2 ■ 

8 Jo JA 

Now, it is enough to take the supremum over G g Cq' 2 ([0,T] x [—1, 1]) to obtain 

< 2^(tt) + ^ j T dt ( (V7r t ) 2 . 

o JO JA 

The inequality in the left hand side of the statement is obtained in the same 
way. □ 
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Setting j3 = in the boundary value problem (2.9) one gets the following bound- 
ary value problem for the heat equation: 

'd t p = Ap in Ax(0,T), 
Po(-) = 7(0 ^ A, (5.8) 
_p t (±l)| = p ± for 0<t<T. 

Lemma 5.8. Let p^ be the solution of (5.8), we have 



W 0) l7) < ^f Q dt(a{p?\j~*{Vp^Y) 



(5.9) 



16 Jo 

Proof. For any G in C 1,2 ([0, T] x [-1, 1]), see (2.11), we have 

J§(P {0} ) = -^£dt(*(pW)V(J™™* P l 0) ) , VG t 
1 " T 



< 
~ 8 



dt(a{pf\(VG t ) 2 ) (5.10) 
f 1 jt/,/„(»)i fY7^ r»~, J°)M 2 \ 



d f (a( P r),[V(J noum ^r)]0 ; 



by inequality (3.12), taking o = 1. The solution of (5.8) belongs to L 2 ([0, T], _ff 1 (A)) 
and its time derivative belongs to L 2 ([0, T], iJ _1 (A)). Therefore, 

/£(Ay) = SU P f? [ dt(a(p° t )V(J™™* P ° t ), VG t 



G€C^([0,T]x[-l,l]) 



\ [ dt(a{p%{VG t ) 2 )} . 



By inequality (3.12), this last expression is bounded by 

^ £ dt(v(p° t ),[V(J™™ *„°)] 2 ). 
We conclude the proof by applying Lemma 3.1 and Jensen inequality. □ 

5.3. 7y(-|7)-Density. In this section we show that any trajectory tt e D([0,T],M), 
with finite rate function, I^ttIj) < oo, can be approximated by a sequence of 
smooth trajectories {ir n : n > 1} such that 

lim tt" = tt in D([0,T],M) and lim 4(7r n |7) = 4(7r| 7 ) . 

Definition 5.9. .A subset A of D([0,T], Ai) is said to be l!^(-\j)-dense if for every 
tt in D([0,T], Ai) such that I^tt^) < oo, there exists a sequence {tt" : n > 1} in 
A such that 7r" converges to tt in D([0,T], Ai) and 7^(7r n | 7 ) converges to I^(tt\^). 

Definition 5.10. Let A\ be the subset of D([0,T], Ai) consisting of trajectories 
tt such that I^(tt\j) < oo and for which there exists 5 > such that tt is a weak 
solution of the equation (5.8) in the time interval [0,6]. 

Lemma 5.11. The set Ai is I^{-\^)-dense. 
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Proof. Fix a path it such that I^(tt\j) < oo and let p^ be the solution of the heat 
equation (5.8). For e > 0, define ir £ as 

p| 0) (-) for0<i<£, 

"*(') = \ /4°-tG) fore<t<2e, 
7r t _ 2E (-) for2e<i<T. 

Since lim e _>o 7r e = 7r in D([0, T], A4) and /(-|7) is lower semicontinuous, it is enough 
to prove that Ve > 0, Z^(7T £ |7) < oo and that lim ,-_ v0 I^(tt £ \^) < I^,(tt\j). From 
Lemma 5.1, for each e > 0, ttq(-) = 7(-) and 7rf(±l) = p±. Decompose the 
rate function X^(7r e |7) as the sum of the contribution on each interval [0,er], [e, 2e] 
and [2e,T]. Since on [0,e] the path ir satisfies equation (5.8), by Lemma 5.8, the 
contribution to the first interval is bounded by 

— F dt ( (X7p[ 0) ) 2 (v)dv. 



This converges to when e I 0. On the time interval [e, 2e], ir £ satisfies 

8t*Z = -dtp^_ t = -Ap&lt = -Arr £ . 
In particular, the contribution to [e, 2e] is equal to 

sup {2 J dt (Vpf \ VG) + I J dt(o-(p { ° } )\7(J ncuTa * p[ 0) ) , VG 

- i [ £ dt(*(pP),(VG t f)}, 



(5.11) 



where the supremum in taken over all G e Go^([0, T] x [-1, 1]). We apply inequality 
(3.12) to the first and second term inside the supremum, then apply Lemma 5.8. 
By Lemma 3.1, the supremum (5.11) is bounded by 



dt I du^y- + / dt(a(p?>)lV(J ncum *pn?) 



^F d tfdu^ + tF dt f^f. 

This last expression converges to zero as e 1 0. Finally, the contribution on [2e, T] 
is bounded by 7^ (7^7). □ 

Definition 5.12. Denote by A2 the subset of A\ of all trajectories 7r such that 
for all < 5 < T. there exists e > such that e < 7r t (u) < 1 — e for (t,u) £ 
[6,T]x [-1,1]. 

Lemma 5.13. 77ie se£ ^2 is 7^(-|7)-dense. 

Proof. By the previous lemma, it is enough to show that each trajectory 7r in Ai 
can be approximated by trajectories in A2- Fix ir in Ai and let p be the solution 
of the equation (5.8). For each < e < 1, let ir £ = (1 — e)ir + ep°. We have that ir £ 
converges to ir as e 1 a.e. in A x (0,T). By the lower semicontinuity of 7^(-|7), 
it is enough to show that 

sup^(7r e |7) < 00 and Urn ^(7^7) < ^(^7) . (5.12) 

e>0 
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Fix e > 0, by construction ttq(-) = 7 and 7if (±1) = p± for almost all t G [0,T]. 
From the convexity of Q(-), for each < e < 1, 

Q(tt £ ) < (1 - e)C(7t) + eQ^ ) < Q(tt) + Q(p°) < 00 . 

Since dt^ e = edtp + (1 — e)<9f7r and, by the assumption, I^(tt\j) is finite, it follows 
from Lemma 5.2, that ir £ G L 2 (0, T; H 1 (A)) and 9 t 7r e G L 2 (0, T; if _1 (A)). Next 
we show that I^ir 6 ]^/) is finite uniformly on e. We decompose the rate i^(-7r e |7) in 
two terms: 

4(7r e |7) < A x +A 2 (5.13) 

where 

Ai=sup{y cft(^ e ,G t ) + y (V7r t e ,VG t ) - i / dt((iK),(VG t ) 2 )} 

(5.14) 

A 2 = su P { | 1 di(a«)V(J ncum *<), VG t ) - iy dt(a«),(VG t ) 2 )} , 

(5.15) 

and the supremum is taken over G G G ' ([0, T] x [—1, 1]). By concavity of cr(-), 
the term A\ is bounded above by 

(l-e)sup{y dt(dm,G t ) + J di(V7r t ,VG t ) - ijf dt(cr(7rt),(VG t ) 2 >} 

f dt(d t p°,G t ) + J T ' dt(Vp° t ,VG t ) - \ J T dt(cr(p t )VG t ,VG t )} . 

(5.16) 

Since p° solves the heat equation, the second line of the last expression is equal to 
zero, while the first line is bounded above by 27j.(tt|7) which is bounded by Lemma 
5.7. By Schwartz inequality, Lemma 3.1 and Jensen inequality 

A 2 <^ £ dt(a(nt)(j nc ™*\V7r!\) 2 ) < £ £ dt jT (|W t | 2 + |Vp?| 2 ) . (5.17) 

By (5.16) and (5.17) we have that sup/y(7r e |7) < 00. 

e>0 

We are now going to prove Hm e _).o Ij-i^ly) — ^r( 7r l7)- By definition of ir £ , we 
have that, for any G G Gq' 2 ([0,T] x [-1, 1]), 

dt(d t n e ,G t ) = (1-e) f dt{(d t ir,G t )+e(d t p°,G t )} . (5.18) 



e sup 



By Lemma 5.5, there exists F G Hq(<t(it)) such that 7r solves the boundary value 
problem (5.1). Taking this into account and (5.18) we can write 

r-T 

dt(d t 7T e ,G t ) 

£ dt{( - V-7r e , VG t ) + (1 - ^(ft) [^V(J ncum * tt) + VF t ] , VG t )} 

y T ^{( - V7r e , VG f ) + (Ff(^),VG t ) + (§<r(7T t e )V(J neum *7r?), VG t }}, 

(5.19) 
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where 



Ff(^) = (1 - e )a(^)[fv(J ncum *^ t ) + V f] - ^K)V(J ncum *<) . 
By the definition of Ij,, see (2.12), we have that 

%(n s h)=su V { / T dt(Ff(^), VG t ) - i f T dt(a(n$), (VG t ) 2 )} =^(^| 7 ), 



(5.20) 

where the supremum is taken over all G <E Cq' 2 ([0,T] x A). The last equality in 
(5.20) holds because <2(7r e ) is bounded for any e > and then (2.15) applies. Since 
for any e > there exists (5(e) > so that <t(tt £ ) > 6(e) we can apply to the first 
term inside the argument of the supremum of (5.20) inequality (3.12) to cancel the 
contribution of the second term inside the argument of the supremum, obtaining 



2 Jo J a vyn) 

On the other hand, from (5.2), 

4(Al) = \J dt{cr(TT t )VF t ,VF t ) . 
Therefore, to conclude the proof, it is enough to show that 



hm - / dt du K = - dt (a(n)VF t ,VF t ) . (5.21) 

By the continuity of a and the definition of Ff (w e ), 

lim 1 ^^ =a(n t )(VF t )\ a.e. in Ax(0,T). 
e^o a(ii £ t ) 

By the convexity of a — > a 2 , the inequality (3.12), the following inequality 

(a + b + c) 3 < 3(a 2 + b 2 + c 2 ), Va, 6, c G K , 
the concavity of a(-) and Lemma 3.1, for any < e < 1, 
(Ff(7r E )) 2 ^ 3^ CT(7r) ^„^ (V7r)2 +(VF)2 ] 



CT(7T e ) 



[a(n) + a( P °)] [ J™ * ((Vtt) 2 + (Vp ) 2 )] 



2 

Therefore (5.21) follows by Lcbcsguc dominated convergence theorem. □ 

Definition 5.14. Denote by A3 the trajectories it € D([0,T],J\4) such that it is 
the solution of the boundary value problem (5.1) for some F € C ' ([0,T] x A). 

The last step is to prove that A3 is Z^(-|7)-dense. We follow the strategy adopted 
in [7]: given a trajectory it in D([0,T], M) with finite rate function J^(7r|7) < 00, 
from Lemma 5.5, there exists a function F in Hq(o(tt)) such that 7r is a weak 
solution to the equation (5.1). Instead of approximating it by a sequence of smooth 
trajectories in D([Q,T],J\4) (cf. [4], [19], [21]), we will approximate F by smooth 
functions (F n ) and we then show that the corresponding smooth solutions (n n ) of 
(5.1) converge to it in D([Q,T],M) and I^(ir n \j) converges to I^,(tt\j). 
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Lemma 5.15. The set A3 is I^{-\^)-dense. 

Proof. In view of the previous lemma, it is enough to show that for each it in A2, we 
can exhibit a sequence {tt„ : n > 0} in A3 which converges to 7r in D([0, T],A4) and 
such that 7y(7r„|7) converges to I^,(tt\j). Fix n £ Ai- Since I^(tt\j) is finite, by 
Lemma 5.5, there exists a function F £ ^(o^tt)) such that 7r is the weak solution 
to the boundary value problem (5.1). We claim hat F £ L 2 ([0, T], H 1 (A)) and 
then, can be approximated by a sequence of smooth functions {F n ) n >\. Indeed, 
let < S < T be such that, it is the solution of the heat equation (5.8) in the time 
interval [0, 5}. We have that VF = -|V(J neum * it) in [0, S] x A and 



f dt f \VF t {u)\ 2 du = [ dt [ ^-|V 
Jo J A Jo J a 4 



(J ncum *7r)(i,u)| 2 du 



T * 

dt / |V^ t (u)| 2 du. 

S J A 



(5.22) 



On the other hand, since ir £ A2, there exists < e < 1 such that e < TTt(-) < 1 — £ 
for 5 < t < T. Therefore 



T 



2 



1 



dt / VFt(u) du < —r-r\\F 



1 2 



^J""" 00 (5.23) 

It follows from (5.22) and (5.23) that F £ L 2 ([0,T},H 1 (A)). Let (F") n>0 be a se- 
quence of functions in G 1,2 ([0, T] x A) such that lim F n = F in L 2 ([0,T],H 1 (A)). 

n— K+oo 

For each integer n > 0, let 7r n be the weak solution of (5.1) with F n in place of 
F. By (5.2) 

= \[ dt {VF? ■ *W)VF?) < dtj^du ||VF t »|| 2 . 

Since the sequence (i 7 ' Tl ) rl >o converges to F in L 2 ([0, T], if 1 (A)) , it follows from 

the last inequality that Ij,(n n \^f) is uniformly bounded. Thus, by Theorem 5.4, 
the sequence ir n is relatively compact in D([Q, T], Ai). Let {7r™ fc : k > 1} be 
a subsequence of 7r' 1 converging to some 7r° in D([Q, T], A4), then {7r" fc : k > 1} 
converges weakly to ^° in Z 2 ([0, T] x [—1, 1]). Since J T (7r n |7) is uniformly bounded, 
by Lemma 5.2 and Lemma 5.3, ir nk converges to ir° strongly in L 2 ([0,T] x [—1, 1]). 
For every G in C ^ 2 ([0,T] x [-1, 1]), we have 



(7r^,G T >-<7,Go) = 



T 







f dt (7r t " fc ),AG t ) -p+ f dtVG t (l) + p- [ dtWG t (-l) 
Jo Jo Jo 

f T ( VG t , CT« fc ) [| V ( J noum * ) + ] ) dt . 
Jo 2 
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Letting k — > oo, we obtain 



(7r°,G T )-(7, Go) 



(7T t °,9 t G t ) 



/ dt(TT°,AG t ) -p+ f d*VG t (l) + p- [ dtVGt(-l) 
Jo Jo Jo 

- j\vG u a(^) [| V( J ncum * ^ t ) + VF t ] ) dfc. 



That is 7T° is a weak solution of equation (5.1). Thus, by uniqueness of weak 
solutions of (5.1), 7r° = n. 

To conclude the proof of the lemma it remains to prove that lim„_ i . 0O Ij,(TT n \j) = 
^ri^ll)- The se q uence ( 7r ™)n>o converges to it strongly in L 2 ([0,T] x [—1,1]) and 
the sequence (F n )„>o converges to F in L 2 ([0, T], iJ 1 (A)). Taking into account 
that 7r is bounded and a is Lipschitz, we obtain 

(■T 



lim 4(ir n \j) 



lim J / dt (VF t n • a«)VF t ") 
-- i y (Vfl • o-(Trt)Vfi) = /|(7r|7) ■ 



□ 



5.4. Upper bound. Let Q = Q 2 be the functional defined in (3.29) with 5q = 2. 
For < a < 1, denote by £ Q : D([0, T], M) -> [0, +oo] the following functional 

£ a {ir) = /|(7r| 7 ) + o(l + a)C(7r) . 
The proof of the upper bound relies on the following proposition. 

Proposition 5.16. Let JC be a compact set of D([0, T], M). There exists a positive 
constants C , such that for any < a < 1, 



Km -LlogQ^JC) < 



inf £ a (ir) + aC(T + l) 



Proof. Fix a density profile 6 : [-1, 1] -> (0, 1), a function G in G 1,2 ([0, T] x [-1, 1]) 
and a function H in G c °°([0,T] x A). For a local function * : {0, 1} Z -> M, c> 
and e > 0, denote £^'* c and -E^ c the set of trajectories (%)t e [o,T] defined by 



G,* 



{r?. g D([0,T],5 JV ) : | j V^(t,r, t )dt\ < c} 
{r?. g D([0,T],5iv) : | / W£(t,»fc)dt| < c} , 



E N,c 

where V^f and W§ are denned in (3.18) and (3.26). 
Define * x (r?) = [r?(l) - r?(0)] 2 , let A^ c be the set 

aG,H _ R VG,*i n p VG n R Vff,*i 

By the superexponential estimates stated in Proposition 3.6 and Proposition 3.7, 
it is enough to prove that, for every < a < 1, 



1 



lim lim lim — lo„ . 



/cn A 



G,H 
N,s,c 



< -- 



1 



1 + a tts/c 



inf £ q (tt) + aC(T + 1) 
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For H e C c °°([0,T] x A), recall from (3.28) the definition of Qh(tt) = Q 2 h (tt), with 
Sq = 2 and write 
1 



— log Q N 
1 



)CnA 



G.H 



logE 



p,N 



By Holder inequality the right hand side of the last equality is bounded above by 



(IT^ 10 ^ 



^logE 

(l + a)N S r > 



l{ICnA^ c }e- a ^ N ^^ 



D (l+a)NQ H (ir N *L e ) 



(5.24) 



From Lemma 3.8, the second term of this inequality is bounded by aC\(T + 1). 
Consider the exponential martingale M t G defined by 



Mi 



exp 



{n[{4 



. G t ) - (tt , G / 



1 

iV 



,G.) {ds + N 2 Cn) e N(^. Gs ) ds j | 



(5.25) 

Since the sequence {rj N : N > 1} is associated to 7, an elementary computation 
shows that on the set A 



G,H 

N,e,c 



exp 



N{ j£(TT N * ie) + Og{e) + 0(C)} , 



(5.26) 



where Og(e) (resp. O(c)) is a quantity which vanishes as e ! (resp. c i 0) and 

Jq{-) is the functional defined in (2.11). Consider the first term of (5.24) and 
rewrite it as 

1 



logE^ [M^(M^)- 1 t{K n A G N ^ c }e 



■a{l+a)NQ H {^ N *i.e) 



(l + a)N 

Optimizing over tt n in JC, since M t is a mean one positive martingale, the previous 
expression is bounded above by 

-— |— inf \jg{w*i e ) + a{l + a)Q H {w * L e )\ + G (e) + 0(c). 
1 + a ireK I ^ J 

Optimize the previous expression with respect to G and H. Since the set K, is 
compact and Jq(- * t e ) and Q//(- * t e ) are lower semi-continuous for every G, -ff, 
e, we may apply the arguments presented in [24, Lemma 11.3] to exchange the 
supremum with the infimum. In this way we obtain that the last expression is 
bounded above by 

- ^ sup inf \-ji(ir*t e ) - a{l + a)Q H (n*i e )} + G (e) + 0(c). 

Letting first e \. and then c \, 0, we obtain that the limit of the previous expression 
is bounded above by 

—L- sup inf {_ jPfr) _ a (l + a )Q ff (^)} 



This concludes the proof of the proposition because sup G Jq{^) = X£(7r|7). 



□ 
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Proof of the upper bound. Let JC be a compact set of D([0,T], M). If for all 
tt e /C, Q(tt) = oo then the upper bound is trivially satisfied. Suppose that 
inf^gx; {S(tt)} < oo, from Proposition 5.16, for any < a < 1, 

Jim" 1 log Q*" (AC) < inf £ a (n) + aC(T + 1) 



-i- inf {4(7T| 7 ) + a(l + a)Q(7r)} + aC(T + l) 



< -- 



1 



— - inf {4(7T| 7 )} - a inf Q(tt) + aC(T + 1) . 
1 + a TreK tts/c 



To conclude the proof of the upper bound for compact sets, it remains to let a J, 0. 

To pass from compact sets to closed sets, we have to obtain exponential tightness 
for the sequence Q^'^ ■ The proof presented in [15, Section 10.4.] is easily adapted 
to our context. 



5.5. Lower bound. In this section we establish the large deviations lower bound. 

The strategy of the proof of the lower bound consists of two steps. We first prove 
that for each ir £ A3, recall its definition in 5.14, and each neighborhood A/^ of it 
in D([0,T],M) 

M llogQ5f{A/;}>-4(7r|7). (5.27) 

The proof of the lower bound is then accomplished by showing, sec Subsection 5.3, 
that for any tt £ D([0,T], M.) with Ij,(n\'y) < 00 we can find a sequence of ir k £ A3 

such that lim fe ^ 00 7r fe = n in D([Q,T],M) and lim fe ^oo It^Ii) = It Ml)- 

The proof of (5.27) is similar to the one in the periodic case, see [15, Section 
10.5.]. It depends on establishing laws of large numbers, in hydrodynamic scaling, 
for weak perturbations of the original process, and controlling by the Girsanov for- 
mula the relative entropies of the processes that go with these perturbations. Fix a 
path ir £ A3- Then, by construction, there exists F £ Cq' 2 ([0,T] x [—1, 1]) so that 
7r is the weak solution of the equation (5.1). Recall from (5.25) the definition of the 
exponential martingale Mf, . Let P^ N be the probability measure on the path space 

D([0,T],S N ) with density M% with respect to P^f: F^ N [A] = E^'f [Af£ 1{A}]. 
Let {f]t)te[o,T] be the process with law P^ N on D(\Q,T],Sn)- Let (7r^) te [o,T] be 
the corresponding empirical measure. Then (^)te[o.T] converges weakly in prob- 
ability to (7r*)tG[o,T]- From the super exponential estimates, Proposition 3.6 and 
Proposition 3.7, we have 



— logQ p ' Ns Kr 1 ^ - ,? ™ — K 

F ItojS.JV 



where H (P^ N | P^w ) stands for the relative entropy given by 



TT 

n \~ v" r ij 



r uir N 



n' 



(II 

To conclude the proof, it remains to show that 
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The proof of [15, Theorem 10.5.4] is easily adapted to our model. □ 

6. Appendix 

In this section we summarize the properties of the equation (2.9) needed to prove 
the main results of the paper. The proofs of these results are based on applying 
standard tools in partial differential equations, although some care need to be taken 
because of the presence of the nonlocal term. Notice that because of the nonlocal 
term the comparison property does not hold for this equation, so tools based on 
maximum principle will not work for (2.9). 

We recall the notion of weak solution of (2.9). A function p(-,-) : [0, T] x 
A — > [0, 1] is a weak solution of the initial-boundary value problem (2.9) if p G 
L 2 ([0,T],H 1 (A)) and for every G G C^QO,! 1 ] x [-1,1]) one has l%{p,i) = 0, 
where £q was defined in (2.8). 

Theorem 6.1. For any ft > there exists an unique weak solution of (2.9). 

The existence of a weak solution of (2.9) is a consequence of the the tightness of 
{Qu n )n>i and the characterization of the support of its limit points, see Lemma 
4.3. The uniqueness can be easily proven performing estimates as in Theorem 6.2 
for all ft. A proof of existence without invoking the hydrodynamic limit can be 
done applying in our setting the argument done in [9], Section 4. 

Theorem 6.2. There exists fto depending on J neum and A, so that for ft < ft 
there exists an unique weak stationary solution p of (2.10). Further, let Pt(po) be 
the weak solution of (2.9) with initial datum po G M.. For ft < fto, there exists 
c(ft) > so that 

\\pt(po) - P~\\lz(A) < e" c(/3)t ||/9o - p\\l*(A)- 

Proof. Let pi t o G A4 and p^t be the solution of (2.9) for t > , with initial datum 
pi o, ? = 1,2. Set v = pi — p2, we have 



1 d 
2dt 



\\vt\\h=- f Vv[Vv-ft X (pi)V(J ncunl *Pi)+ftx(P2)V(J ncunl *P2)] 

J A 

= - f Vt-[V«-/?x(pi)V(J ncum * U )+/?[ X (p 2 )-x(pi)]V(J ncum * /92 )]. 

J A 

(6.1) 

Since X (a) < j for a G [0,1] and \ X (p 2 ) - x(pi)\ - \{p2~pi} [1 - (p2 + pi)}\ < \v\ 
we have that 

j Vv [Vv - ftx(pi)^7(J ncum * «)] > ||V«||i a [l - |] (6.2) 
/ V« [x(p 2 )-x(pi)]V(J ncum * P2 )] 

J A 



<sup|VJ noum |||V7;|| La ||t;|| La 



<|sup|VJ" cum |1|V V ||i 2 + ^|k||i 2 , 



(6.3) 



for any a > 0. Taking into account (6.2) and (6.3) we can estimate (6.1) as following: 
~\Mh < -||V«||| 2 [1 -ft{\ + a - sup |VJ—f)] + A|M| 2 L2 . (6.4) 
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and we choose a so that 



i + ^sup|VJ"°H 2 <i 



Since we are in a bounded domain we can use the Poincare inequality 

\Hh < c{k)\\vv\\l- 2 

obtaining 
Take (3q so that 

[i-f] ft n 

(7(A) 2a 

Then for /3 < f3 there exists c(/3) > 0, - ^ = c{fi) > 0, so that 

~M\h<-m\Hh' (6-6) 

This implies immediately that the stationary solution is unique and that it is ex- 
ponential attractive in L 2 . □ 
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